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V. A. Sandborn 


SUMMARY 


The problem of dust in the special environment of Spacelab 
is reviewed, A major factor of importance in the zero-gravity 
environment will be the presence of large particles. It is neces*- 
sary to consider aerosol dynamics at Reynolds numbers larger than 
those of importance in the earth environment. The drag coefficient 
of particles in the range of Reynolds number from 1 to 100 is 
examined in detail. The accurate relation between the drag coef- 


ficient, c p, 
found to be 


and Reynolds number, R , for the range I^R-^IOO was 

S Q 



Based on drag data for spheres currently available for slip and 
free molecular flow a general graphic correlation for the drag 
coefficient as a function of Reynolds and Knudsen number was 
developed. 

The general equations that govern the dust dynamics were 
reviewed. 


I 
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INTRODUCTION 

Particle dynamics associated with the special environment 
of Spacelab may pose a new set of problems . The near zero 
value of gravitational forces will greatly reduce sedimenta- 
tion of both large and small particles. The problem is further 
augmented by the presence of "g-jitter" effects associated 
with the space craft. Large particles, which would be created 
in an industrial use of the zero-g environment, could pose a 
major difficulty. Under the influence of gravity the large 
"dust" particles created, in for example a machining process, 
quickly settle out and do not pose a problem. In the Spacelab 
environment the large particles would float free and would also 
detach from solid surfaces easily. The slight variation in 
gravity (g-jitter) associated with Spacelab would act to keep 
the large particles in suspension. 

A great deal of information, both experimental and theoreti- 
cal has been obtained on the dynamics of fluid-solid suspensions, 
ref, 1 through 5. Attention has been focused heavily on the dy- 
namics of small aerosol particles, since the newly developed laser 
velociiheters employ aerosols as tracers to measure fluid motion. 
Problems associated with aerosol particle motion in turbulent 
fluid flows were considered in detail by Lumley, ref. 6. In 
general the small aerosol particles can follow approximately the 
local turbulent motion, so that detailed information on the tur- 
bulent motion is obtained with the laser systems. The small 
particles would not pose a major problem in Spacelab, as they 
could be trapped on filters in a common air cleaning process. 

The large particles, because of their inertia, will not respond 
as quickly to the local fluid motion. Thus, it may not be pos- 
sible to filter out the large particles as efficiently as the 
small particles in the Spacelab atmosphere. Either a larger air 
flushing system in Spacelab, or limitation on dust generation, 
may be required . 
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GENERAL IDENTIFICATION OF AEROSOL PARTICLES 

Particles of interest in the* spacelab environment will most 
likely be generated by the disintegration of solids. There is 
also the possibility of gas-to-particle generation in industrial 
application. Liquid particles might also be generated in cer- 
tain experiments. Each of these types of particles will be 
referred to as aerosols. The aerosols will be assumed to be 
inert with respect to the spacelab atmosphere. In general, 

Hidy and Brock, ref. 3 $ have pointed out that the first criter- 
ion for defining an "aerocolloidal system 1 ' involves a "small" 
rate of settling of the suspended particles. Thus, the present 
discussion will fall in the broad spectrum of aerocolloidal systems. 

Particle Size - The jargon employed in aerosol science refers to; 
dust as solid particles produced by disintegration , smoke or fumes 
as particles formed from gases, and mist as liquid droplets. A 
very wide range of sizes and masses exist for aerosols. Friedlander, 
ref. 3> estimates that aerosol sizes can vary by factors of 10^, 
while the mass can vary by factors of 10^. For the special environ- 
ment of spacelab the range of sizes and masses may be even greater. 
The minimum size of an aerosol will be that of a "molecular cluster", 
which is of the order of 10$. Obviously, the very small molecular 
clusters are not likely to be as important in the present discussion 

as the larger dust particles. Particles formed directly from the 

—6 

gas state are usually less than a micron, in size (1 micron =10 
meters, and is usually written as 1jw.no). Dust particles range from 
Q.l/uwup to 1 00/im f or atmospheric type cases. Again the large "dust" 
particles generated in machining processes are greater than lOOQju-m 
in size. 

If the particle is spherical in shape it is obvious that the size 
or characteristic length will be the diameter. If the particle is 
not spherical then the characteristic length will be related to a 
aerodynamic dimension , The aerodynamic dimension, or equivalent 


4 


"effective” diameter would be defined based on the drag of the 
particle-* The effective aerodynamic diameter is determined by 
use of a cascade impactor. The cascade impactor is one instru- 
ment commonly used to classify aerosol sizes. The instrument con- 
sists of a series of collector plates, as shown on figure 1. Both 



Figure 1 . Cascade impactor. 

the sizes of the orifices and the clearance between the orifice 
and the collector plates are varied from one stage to the next. 

The gas flow velocity is progressively increased for each stage 
so that smaller and smaller particles fail to turn with the flow 
and are collected on the plates. Any number of stages may be em- 
ployed to obtain a measure of the size distribution of aerosols . 
The efficiency of a stage is a function of the Stokes number 

c >o Ud 2 
S = ? P 

t 1&u.d (1) 

where 

c slip correction factor 

s 

S o ^ particle density 
U gas velocity through the orifice 

d is the effective particle diameter 

P 

M. is the gas viscosity 
d is the orifice diameter 

Thus, the effective aerodynamic diameter can be defined as the 
diameter of a hypothetical sphere of unit density with the same 
Stokes number as the particle in question. The efficiency of the 
stages of a cascade impactor is such that a range of particles 
will be caught. Thus, in practice the stage is usually identified 
by a Stokes number which corresponds to the 50% impation efficiency 
level . 
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Cadle, ref. 7, describes a number of different methods that 
have been developed to define particle sizes. One technique that 
is the "equivalent ” circle whose area is equal to the projected 
area of the particle. Special graticules have been developed to 
allow the particles to be evaluated under a microscope. For the 
present discussion it would appear that the aerodynamic specifi- 
cations of particle characteristic size will be the logical para- 
meter. 

The aerosols of interest will be polydispersed in size with a 
range far greater than that encountered in a gravitional environ- 
ment. The physical description of importance will be the size 
distribution function. For a given position in space the number 
of particles per unit volume can be determined if the size dis- 
tribution function, n, is known. The number of particles, N, 
with' sizes in the volume range V to V+dv is 

AN = ndv (2) 

For near spherical particles the distribution function might better 
be expressed in terms of the effective diameter, dp rather than the 
volume. The relation between the volume and diameter distribution 
function, n , is, (Friedlander , ref. 5) 

n -TC^p 11 (3) 

n<i ~ 

Development of the time history of the particle distribution func- 
tion within spacelab ?/ill be necessary to define the specific dust 
problems . 

Four specific effects can be identified in the time history of 
the distribution function: 

1. IT xtcleation - formation of particles from the gas phase. 

2. Coagulation - clumping of particles. 

3. Sedimentation - fall out of particles due to gravitational 
forces . 

4. Diffusion - 3rowian motion and flow convection produces 
diffusion of the particles. 
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Nucleation and coagulation are processes which occur within the 
volume, dv , of gas, while sedimentation and diffusion represent 
transport across the volume boundaries. For the application to 
spacelab, sedimentation must be viewed as a highly time dependent 
phenomenon, due to the "g-jitter" aspect of the craft. Within 
the global aspects of the spacelab problem, the "sticking’ 1 of 
particles to the surfaces must also be considered. While small 
particles may adhere to surfaces they come in contact with, the 
large particles will easily pull away when changes in acceleration 
occur . 

Friedlander, ref. 5, suggests the distribution function be 
divided as follows 


dn _ [dja] f— 1 

dt “ LdtJi + LdtJ 


(4) 


where i corresponds to the processes internal to the volume and e 
corresponds to the processes external to the volume. The moments 
of the distribution function 


N 


\ = /“X d < d n> 


n J d p ' p' 

can be related to specific physical properties: 


(5) 


1 


4 . 


7 . 


M gives the total concentration of particles in suspension 
at a given location. 

divided by M q is the average diameter, d^ f of the particles 
in suspension. 

ML times 7X gives the total surface area per unit volume of 
particles, and the average surface area is given by 11 ^2/M 

o 

l-U times 11/6 gives the total volume of the particles, and 
the average volume is obtained by dividing by M Q . 

is proportional to the total surface area of the particles 
that sediment from the volume of fluid. 

Me- is proportional to the mass flux of particle sedimentation 
from the fluid. 

is proportional to the scattering of light from particles 
that are smaller than the wave length of the incident light . 
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Obviously the larger particles will contribute greater, or dom- 
inate the higher moments of the distribution function. 

The possible shapes of the size distributions are reviewed by 
Cadle, ref. 7. A common occurring probability function for dust 
particles in nature is the log-normal distribution. The log- 
normal distribution is skewed, in that the probability of finding 
very small particles drops off rapidly, while the probability of 
large particles occurring is more likely. The skewed distribution 
is expected to occur in spacelab since the lack of sedimentation 
will act to increase the number of large particles in suspension. 
Also the very small particles will be more likely to stick to ex- 
posed surfaces and thus diffuse out of the atmosphere of spacelab. 

DEVELOPMENT OF THE EQUATIONS GOVERNING DUST DYNAMICS 

A general relation for the time rate of change of the size dis- 
tribution function, n, can be developed from equation (4)* Formu- 
lation, coagulation and grov/th are processes internal to a given 
volume, while diffusion and sedimentation involve transports across 
the walls of the volume. 



Figure 2. Elemental Volume. 

External Processes - For the volume a, b, c, d, e, f, g, h, 
shown In figure 2 the net accumulation of particles within the 
element due to convection Is 

Convective Accumulation = -<5x6ySz|T“ + (6) 

where u, v, and v are the components of the convective velocity 
vector V. The accumulation of particles due to diffusion is 

Diffusion Accumulation = Sx6y5zD^~^ + (7) 

where D is the diffusion coefficient. Equation (?) is for an 
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isotropic fluid with equal diffusion in all three directions. 

Sedimentation is related to an external gravitational force 
field, which will be very small and varable for the spacelab , 
environment. Other force fields, either electrical or thermal, 
may also be present in spacelab of equal strength to that of 
gravity. If a force field is present the particles will migrate 
in the direction of the force with a velocity, c, given by 

c =£ ( 8 ) 

f 

where F is the force and f is the resistance to motion of the 
particles due to fluid drag or friction. The accumulation of 
particles due to a force field can be written as 

Force Field Accumulation =*6x6ySz(c n+cn+cn) (9) 

x y z 

The rate of particle accumulation in the volume 6 y cz.. due 

to convection, diffusion and external force fields can be written 
as 

(|“) e = -r-STn + DV^N -V?h (10) 

In this equation the velocity, v, is for the most part indepen- 
dent of the particle concentration ■, and dependent on the fluid 
mechanics of the system considered. For spacelab application 
the velocity would be due to the air exchange system within the 
working area. Obviously, if large particles must be drawn into 
the air cleaning filters, then the magnitude of T will be large. 
The flow within the working area will be complex and not easily 
predicted, due to the obstructions and personnel. 

The diffusion coefficient , D, and the drift velocity, C, will 
both depend on the particle size. The diffusion coefficient de- 
creases as the particles become large. For a constant force 
field, which is independent of particle size, the drift velocity 
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will decrease as the particles become larger. The friction fac- 
tor, f, in equation (8) is directly proportional to the particle 
diameter for the small particles. In the "Stokes" flow region 
the friction factor for a sphere is 


t ~ 3 W d p 


( 10 ) 


where m is the fluid coefficient of viscosity, d^ is the particle 
diameter and C$ is a slip correction factor. For large particles 
above the Stokes flow regine., the friction factor may be less 
dependent on diameter than that given by the Stokes relation. 

Obviously there are many cases when the force field will also 
be a function of the particle size. Gravitational forces will 
increase as the mass of the particle. Thus, it is difficult to 
specify how the drift velocity will vary with particle size. 

Internal Process - The internal processes of growth and coagu- 
lation of particles have been treated in detail by Friedlander, 
ref. For growth the net rate of change in particle number 

in the elemental volume €v — Sx Sy <fz is, ref. 5, 


a i 

*t ; g "“av 


OD 


Which is a particle continuity relation. I(v,t) is the particle 
flow "current" or the number of particles per unit time per unit 
volume of gas passing a point Sv . 

ifflen particles collide some will adhere together resulting in 
larger particles. Coagulation produces large particles from 
the smaller particles, which results in the skewed log normal 
type of particle distributions. Friedlander, ref. gives the 
following formulation for the coagulation 

( 12 ) 

(■“TV = Formation - Loss 

dt coag 

When any two particles of volumes \J , and adhere to form a new 
particle the number density of particles is reduced by one. 
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Obviously, the formation rate will depend directly on how often 
particles collide. Employing a "collision frequency function", 
jV x ) such that the collision rate can be expressed as 

Collision Rate -/a (v, } V z ) n(v ; ) r\(v 2 ) d V, dv 2 
Thus the formation of particle of size V, , from the collision 
of smaller particles of size ( V, - V z ) and is 

Formation over dv, = ^ v,-v» ) »<VJ nfv; -vjdvj 4 v, 

and the loss of particles of size V, , due to collisions with 
other particles is 

Loss Over dV^ - fj*# V z ) n CV,) V>( V x ) d VJ cl ’ 

The net rate of formation of particles given by equation (12) is 


y 

(f£) = if ^(v 2 ,v,-v z ) Tn(vjri(v,-v Odv* 03) 

-C#(v tJ vO n(v x )^(v,)dv^ 

Dynamic Equation «■ The general dynamic equation for the size 
distribution function is 

y 

+ V *70 “ 7 - D Vr) + zf 0 /&(v 9 j V,- Vy) Tn(v z )r)(v,-Vi)<lv z 

-fs8(v li v^)ri(v,)Tn(v 9 )<iv*--7'Cy) (1Zf:> 


For the case of incompressible flow, which would be the flow en- 
countered on the space laboratory filtering system 


V - m — V • V ~n 


(15) 


AEROSOL DRAG FORCES 

Removal of dust from an environment depends mainly on the 
movement of air through a filter system. For the specific space- 
lab environment filters will be of major importance, as gravita- 
tional sedimentation will not be present. Thus, the drag forces 
on the aerosols will be the basic mechanism for the management 
of dust in the spacelab. The general view of external processes 
considered in the previous section employed a friction factor, f, 
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to represent the drag force on the particles . It was also noted 
that the size distribution of particles can be expressed in terras 
of a Stokes number, eq. (1), which was directly related to the 
particle aerodynamic drag. In this section a review of the drag 
process is covered, and the areas of specific interest to space- 
lab are identified. 

The size and shape of dust particles that will be present 
in the spacelab environment can not be specified. As noted 
earlier, it is obvious that much larger sized particles can be 
expected due to the absence of sedimentation. It is also pos- 
sible that some reduction in the number of very small particles 
will occur due to their attachment to the very large particles- 
not normally present in gravitational force fields. Although 
the aerosols found in spacelab are not expected to have a well 
defined shape, the present discussion will consider first the 
ideal spherical particle. As noted in the discussion on size 
distributions it will be convenient to describe random shaped 
particles in terms of an equivalent "aero dynamic M spherical dia- 
meter . 


Sphere Drag - For the discussion of drag it is convenient to 


employ a non-dimensional drag coefficient defined as 


u (n Ub; 

where F(=cf) is the force on the sphere, as employed in equation 
(8 ), yo is the air density, 1/ is the air velocity, and 7TA Z is 
the frontal area of the sphere. The drag coefficient is found 
to be a function of the non-dimensional parameters, Reynolds 
number, Knuds en number and in the general case the Mach number 
also. The Reynolds number expresses the ratio of the inertia 
force to the viscous force of the flow. For the spacelab appli- 
cation the Reynolds number will be the main factor in the deter- 
mination of the drag on the aerosol particles. The Knuds en 
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number is the ratio of the mean free path between molecules of 
air to the sphere diameter. For large Knudsen number the drag 
approaches the free molecular flow value, while small Knudsen 
number flows produce continuum flow results. Aerosol problems 
normally are found to contain many particles in the region be- 
tween continuous and free ;molecular flow. For spacelab consider- 
ations the possibility of a great number of large particles may 
make Knudsen number considerations of only limited interest. 

The Mach number is the ratio of the flow velocity to the speed 
$>f sound. Mach number effects are limited to high velocity flows 
in the continuum regine . Mach number effects are not important in 
the spacelab dust management considerations. 

Figure 3 shows the values of both Knudsen and Reynolds number 
as a function of sphere diameter. The environment of spacelab 
will be nearly that of a standard atmosphere, with a mean free 
path of 0.063 microns. Reduction of the pressure or more 
likely an increase in temperature will produce slight increases 
in the mean path. The curve labled 1/2 atmosphere represents an 
upper limit on the variation of Knudsen number that might be ex- 
pected in the spacelab environment. Knudsen numbers of the order 

of 3 to 10 are required for free molecular flow. In the range of 

-3 

Knudsen numbers between 3 and 10 the flow is referred to as 
"slip flow" , since the flow velocity is not zero at the sphere 

_3 

surface. For Knudsen number less than roughly 10 the flow can 
be treated as a continuum, and no slip occurs at the sphere surface. 

The Reynolds number determines whether the flow around the 
sphere is dominated by the viscous forces or the inertia forces. 

For Reynolds numbers less than roughly 0.1, inertia forces can be 
neglected completely. The classical Stokes solution for the 
drag of spheres, eq. (10), is for the case where inertia terms 
are neglected. The Stokes solution in terms of the drag coef- 
ficient and the Reynolds number may be written as 



( 17 ) 



Knudsen Number 
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Figure 3. Variation of Knudsen and Reynolds numbers with sphere diameter. 

For Reynolds numbers below 0.1 the relation is very close to the 
actual measured drag coefficients. The relation is roughly 1 0 % 
too low for a Reynolds number of 1 • Oseen, ref. 8, employed a 
simplified model for inertia to obtain a second order approxi- 
mation to the sphere drag at low Reynolds numbers 

(,8) 

D Re 1b 

This solution improves the prediction of drag up to a Reynolds 
number of roughly 0.5. It over predicts the drag by 7 % at a Rey- 
nolds number of 1 . Further development of the Oseen approach 
lead to the inclusion of higher order Reynolds number terms, how- 
ever the improvement over the predictions of equation (18) is very 

At Reynolds numbers greater than 1 the viscous boundary layer 


minor . 


Reynolds Number 


developed on the sphere starts to alter the flow around the 
particle. Above a Reynolds number of approximately 40, separa- 
tion of flow from the rear surface of the sphere becomes pro- 
nounced. At much higher Reynolds numbers , outside the area of 
interest for aerosol problems, the viscous constraints become 
small compared to the inertia terms and turbulent flow developes. 
In the area of interest to the spacelab problem, the Reynolds 
numbers are not expected to be greater than 100. 

5* or the drag of spheres in the range of 1<R <400 the empiri- 
cal relation proposed by Klyachko, ref. 9* is found to fit the 
experimental results within 2% 

c/ 3 

(19) 


2/ 3 

c +22_3 

C D ~ Re' 1 + 8 3 


For the spacelab study it might appear that the interest will 
be for Reynolds numbers less than 100, fig. 3. A comparison 
of equation (19) with measured values of (compiled by Fucks, 
ref. 2) is shown in figure 4. In the Reynolds number range from 
approximately 3 to 100, equation (19) is slightly low. A some- 
what improved fit to the experimental can be obtained by adjust- 
ing the constant of equation (19) to 

0.7 


c _ + 2® ) 

D ~ Re' 1 6.8 ; 


( 20 ) 


Equation (20) is also plotted on figure 4. Equation (20) gives 
a slightly more accurate comparison with the data in the range 
of interest for spacelab conditions. 

As noted previously the drag depends not only on the Reynolds 
number, but also on the Knudsen number. The data of figure 4 
is for the continuum flow conditions where the knudsen number is 
very small. As noted in equation (10), if the Knudsen number is 
of the order of 10 1 or greater it is necessary to employ a slip 
flow correction for the drag coefficient. At this point it is 
not obvious that the small particles, which correspond to the 




Drag Coefficient 




Number 
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larger Knudsen numbers, will be of major importance in space- 
lab conditions. 

The "slip" flow correction for the drag of a sphere, C c , was 
developed by Millikan, ref, 10, 

r_ 

= /+£/<* (fl-t-Be'**") (21 ) 

where A = 1 .257, B = 0.2*00 and C = 0.55, ref. 11 . 

Figure 5 was developed, using equation (19) divided by equation 
(21), to determine the different Knudsen number curves. For 
the spacelab application the slip flow conditions are probably 
limited to Reynolds numbers less than 10*^ , However, to define 
the expected extent of the complete sphere drag map the limiting 
case of large Kach number drag measurements, ref. 12 and 13, is 
also plotted on figure 5. For large Mach numbers (or more im- 
portant to free molecular flowj large speed ratios) a limiting 
value of is applicable for Knudsen numbers greater than 5 
to 10. For smaller Knudsen numbers it is not possible to 
theoretically predict the limiting value of C^. Experience 
with heat transfer measurements for spheres, ref. 14, and cylin- 
ders, ref. 15, suggest that the results for Mach numbers greater 
than roughly 2 gives a limiting lower value for C^. No doubt 
the Knudsen number curves may fair more smoothly into the high 
Mach number measurement, although the region of intersection 
between the curves drawn are at the correctly measured Knudsen 
number . 

Also shown on figure 5 are the Stokes, eq. (17), and Oseen, 
eq. ( 1 8 ) , relations for the low Reynolds number continuum flow. 

As pointed out by Fucks, ref. 2, the Oseen correction to the Stokes 
relation is only usable in a very limited Reynolds number ratvgse. 
The emprical relation of Klyachko, ref. 9, or the relation given 
as equation (20) can better approximate the drag at ’•high" Rey- 
nolds number than either the Stokes or Oseen relations . 



The free molecular flow limit shown on figure 5 corresponds 
to the upper limit of speed ratios obtained from theoretical 
considerations, ref. 16. This limit will not be encountered 
in the environment of space lab . 

Non-Spherical Particles. - The spherical particle is at 
best a highly idealized model of the aerosols that will be 
encountered in the space lab .environment . It is not possible 
to develope a drag "map” for each shape of particle encountered. 
Only the sphere and ellipsoids have been evaluated in detail. 

The ellipsoid can include the limiting shapes of a. cylinder and 
at the other limit a thin elliptical plate or disk. Studies 
of ellipsoidal particles have received considerable attention, 
ref. 2. 

Somewhat surprising it is found that the elliptical particle 
in a gravitational field will not have a preferred direction 
of orientation in a low velocity, uniform, viscous flow. This 
result, that the torque acting on an ellipsoid of revolution 
will be zero, has been observed both experimentally and analyti- 
cally. Obviously particles which are not symetric in shape 
will develope torques which will result in a preferred direction 
of orientation. In a gravitational free fall, elongated particles 
orient themselves in a way that balances the troque due to drag 
with the gravitational force. As a result, the orientation that 
produces either a minumum drag or torque may not be the preferred 
orientation in a gravitational free fall. In the absence of 
gravity, as encountered on spacelab, the orientation will depend 
more heavily on the torque developed by the drag forces. The gf- 
Jitter encountered in the gravitational field of spacelab will 
produce a random force on the particles, which could cause an 
oscillation in the orientation. In the absence of a constant 
gravitational force erratic motion, such as spiral, zigzag or 
gliding trajectories of the aerosol particles would not be en- 
countered . 
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The theoretical evaluation of the drag of ellipsoids for the 
purely viscous flow regine results in a solution similar to the 
Stokers solution for a sphere, eq. (10). Only the numerical 
coefficient is changed to account for the change in shape, ref. 

17. A numerical coefficient, D , which is a function of the 

s 

shape of the particular spheroid, is used to correct the Stokes value. 
For a prolate ellipsoid (rotation of an ellipse about its 
major axis) with motion along the polar axis, the shape co- 


efficient is 


D s- 


4-rdf- jj 


(22) 


where a/b is the ratio of the major to the minor axis of the 
ellipsoid. For the probate ellipsoid with motion transverse to 
the polar axis the shape coefficient is 


»s = 




(23) 


J +(%) 

[(%) " / J 9 - 

For an oblate ellipsoid (rotation of an ellipse about its minor 

axis) with motion along the polar axis the shape coefficient is 

fa, \2 






^)a a / b f-7] , ^ + 

For the oblate effipsoid with motion transverse to the polar axis 
the shape coefficient is 

D _ j-[fe) - -0 

Figure 5 is a plot of the drag coefficient for a number of the 
possible ellipsoid shapes computed from equations (22) through 
(25) . The equatorial axis of the ellipsoid is the characteristic 
used to define the Reynolds number. The oblate ellipsoid pro- 
duces drags less than the sphere, while the prolate ellipsoid 
produces a drag approaching that of a disk. 
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Figure 6. (Concluded) Drag coefficient of ellipsoid shapes 
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For orientation of the ellipsoid at angles rather than parar* 

llel or transverse to the polar axis the drag can be resolved 
into components. If the particle is rotating (due to Brownian 
motion) it can be shown that the average resistance is equiva- 
lent to the polar axis being oriented parallel to the flow one- 
third of the time and perpendicular to the flow two-thirds of 
the time, ref. 2. A statistical mean resistance can be computed 
using the 1 /3 and 2/3 rule . 

In the discussion of particle size it was noted' that an effec- 
tive diameter, d , could be defined such that the drag of the 

particle corresponds to that of an equivalent sphere. For the 

a 1 /3 

ellipsoids the equivalent diameter is d (•£•) J for the prolate 

a —1 /3 e k 

and d Q (~) for the oblate case, where d g is the equatorial 

diameter. Fucks, ref. 2, suggests the use of a "dynamic" shape 
factor, which is the square of the ratio of the equivalent di- 
ameter of the particle to that of a sphere of the same volume. 

The dynamic shape factor is similar to the use of a Stokes number, 
eq. (1), to define particle characteristics. Listed below are 
typical dynamic shape factors for a number of characteristic 
particles . 

TABLE 1 . TYPICAL DYNAMIC SHAPE FACTORS 
(FOR A RATIO OF THE HEIGHT TO THE BASE; DIAMETER, OR RATIO OF 
AXES, of if .00) 

CYLINDER - HORIZONTAL 1 .32 

VERTICAL 1 .07 

PARALLELEPIPED WITH SQUARE BASE 
HORIZONTAL 1 .31 

VERTICAL 1 .07 

ELLIPSOIDS OF ROTATION 

HORIZONTAL 1 .28 

VERTICLE 1 .36 

TWO CIRCILAR CONES JOINED AT THEIR BASES 


HORIZONTAL 


1 .27 
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GLASS SPHERES JOINED IN A CHAIN 


TWO SPHERES 1.1.6 

THREE SPHERES 1 .31 

POUR SPHERES 1 .70 

EIGHT SPHERES 2.14 

THREE FLAT DISKS 1 .26 
SEVEN FLAT DISKS 1 .70 
SIX OCTAHEDRA 1 .31 

SPECIFIC BODIES 

OCTAHEDRON 1 .06 

CUBE 1 .07 

TETRAHEDRON 1 . 1 8 


With the exception of a few of the large chain particles it 
appears that the variations in effective diameter is less 
than 20%. Thus, use of the basic sphere drag information with 
slight corrections for particle shape should give accurate re- 
sults . Because of the wide variety in particle shapes and 
orientations it is impossible to suggest a specific statisti- 
cal average for the shape factor. For spacelab applications 
the main interest will be in moving the aerosol particles to 
the filter system. Thus, for conservative calculations the 
drag coefficient should be selected somewhat smaller than 
expected . 

The discussion of drag of aerosol has been limited to the 
case of rigid particles. If the particle is liquid it may de- 
form or develope a circulation which alters the resistance to 
movement through the air. For a spherical drop with circulation 
the Stokes drag, eq. (17) is reduced by a factor, ref. 18 

i + P'VW) 

(I (26) 

where m is the air viscosity and is the liquid droplet vis- 
cosity. For water in air the correction given by equation (26) 
is less than 1 percent. The liquid must be of low density for 
the circulation correction to be of importance. 
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Very large liquid droplets can be expected to breakup into 
smaller droplets that are stable • Obviously the stable drop- 
let size is determined by the surface tension of the particular 
liquid • 



25 


CONCLUDING REMARKS 

The najor problem associated with dust In the Spacelab zero- 
gravity environment will be the failure of large particles to 
settle out of the air* Also the exlstance of g-jitter of the 
laborator will act to further keep the larger particles In sus- 
pension* Thus, It is necessary to examine the particle dynamics 

in the Reynolds number range from 1<R<100* An examination of 

e 

the measured data in this Reynolds number range led to the develop- 
ment of an emprlcal curve fit for the drag coefficient, c Q , of 
the form 



A survey of experimental measurements of the drag of spheres 
over a wide range of flow conditions was evaluated. A graphic 
correlation of the drag coefficient as a function of both Reynolds 
and Knudsen number was developed* It was possible to cover a 
range of Knudsen number from zero to approximately 50* The pre- 

sent analysis was: Halted to a Reynolds number range between 10 
2 

and 10 , Although information was lacking on the correlation of 
drag coefficient with Mach number at subsonic conditions, It was 
possible to Include the limiting case of large supersonic Mach 
numbers « 

The general equations governing the dust dynamics were reviewed. 


I 
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